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Abstract 

The chaotic properties of the three-site antiferromagnetic Ising model on Husimi 
tree are investigated in magnetic field. Macroscopic quantity of three-site antifer- 
romagnetic Ising model is generated by one dimensional map. It is shown that in 
certain parameter setting strange attractors of this map exhibit multifractal scal- 
ing. By applying thermodynamic formalism we find nonanalyticity in free energy 
as well as in entropy. We show that the temperature of phase transition depends 
on parameter of model and varies from negative to positive value. 
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1. INTRODUCTION 



One of the formalism to analyze nonlinear physics having complicated fractal objects 
and strange attractors is the thermodynamic approach]!], 0]. The advantage of this for- 
malism is that the number of degrees of freedom is usually enormously large and their 
properties are described by thermodynamical function which contain most of the rele- 
vant information about macroscopic systems. In order to characterize these macroscopic 
systems order parameter has been sought. One of this order parameters is Lyapunov ex- 
ponents. The tool with which we can measure Lyapunov exponents and other dynamical 
quantities is similar to the thermodynamic formalism in equilibrium statistical mechanics. 
By using thermodynamic formalism many studies have been done on strange attractors 
and it has been found nonanalitisity in free energy function and scaling in the distribu- 
tion of Lyapunov exponents^, ^. Besides there are many other relationships between 
Ising-like systems and thermodynamic formalism[^]. 

In this paper we investigate the three-site antiferromagnetic Ising spin model (TSAI) 
and show how to obtain connection between the TSAI model and thermodynamic formal- 
ism. The reason for studying the Ising model with multisite interaction is that it plays an 
important role for investigations of real physical systems, such as binary alloys 0, clas- 
sical fluids solid ^He 0, lipid bilayers 0, and rare gases |]10[- Recently the multisite 



interaction Ising model on Husimi tree has been investigated. It is shown, that this ap- 
proach yields qualitatively good approximation for the ferromagnetic phase diagrams than 



conventional mean-field calculation [|IT|. In contradiction to ferromagnetic case, when we 
change the sign of three-site interaction the situation changes drastically. In certain values 
of temperature and magnetic field the TSAI model has nontrivial thermodynamic limit 
and as a consequence the magnetization exhibits chaotic behavior |T2|, |r^. In this case the 
magnetization is not an order parameter. In this paper we introduce Lyapunov exponent 
as the order parameter in chaotic region using thermodynamic formalism. 

With the "thermodynamic formalism" we investigate TSAI model in the chaotic re- 
gion and describe its chaotic properties via the invariants characterizing fractal sets (e.g. 
strange attractors). In particular, we obtain the entropy and free energy function and 
focus on whether the thermodynamical quantities have phase transition 0]. It is, in gen- 
eral, hard to determine such behavior unambiguously by numerical methods if one does 
not have further arguments or exact solutions (as in the case of a; — > 4a;(l — a;)). However, 
there is a possibility to remedies these deficiencies if one considers the characteristic Lya- 
punov exponents as an order parameter, which will differ in two phases. In this paper we 
numerically calculate this order parameter as a function of the "temperature" and show 
where the phase transition occurs. 

The contents of the remainder of this paper is as follows. TSAI system on Husimi 
tree and recursion relation is presented in Sec. 2. In Sec 3. there is a discussion of the 
property of TSAI model in thermodynamic limit. In Sec. 4 we present the TSAI model in 
the case of "fully developed chaos" and obtain the exact connection between this model 
and chaotic attractors. By using the "thermodynamic formalism", the phase transition 
in terms of Lyapunov exponents, free energy and entropy function is analyzed. Finally, 



2 



in Sec. 6 we summarize our results and comment on their implications for the study of 
other systems. 



2. TSAI MODEL ON HUSIMI TREE AND RECURSION RELATION 
The advantage of the Husimi tree is that for models formulated on it, an exact recursion 



relations can be obtained. The pure Husimi tree [Tj], shown in Fig.l, is characterized by 
7, the number of triangles which goes out from each site and n, the number of generation 
. The Oth-generation is a single site from which come out 7 triangles. All subsequent 
generation comes out by gluing up 7 — 1 triangles to each free sites of previous generation. 
The TSAI model in a magnetic field is defined by the Hamiltonian one 

H = -j'^^OiOjOk- fi^ai, (1) 

A i 

where (Xj takes values ±1, the first sum goes over all triangular faces of the Husimi tree 
and the second over all sites. Additionally, we use the notation J3 = /SJg, h = ph , 
(3=1/ kT, where h is the external magnetic field, T is the temperature of the system and 
J3 < corresponding to an antiferromagnetic coupling (in all our numerical calculations 
we put J3 = —1). 

The partition function will be written as 

Z = E ea;p < J3 E CTj ajak + /i ^ (Tj > , (2) 

{^} [ A i } 

where the summation goes over all configurations of the system. 

When the Husimi tree is cut apart at the base site, it separates into 7 identical 
branches. The partition function can be written as follows: 

Zn= Yj exp{hao} [5'„(o-o)]'^~\ (3) 

where cxo are spins of base site, n is the number of generations [n ^ 00 corresponds to 
the thermodynamic limit). Each branch, in turn, can be cut along any site of the 1th- 
generation which is the nearest to the central site. The expression for gn{o'o) can therefore 
be rewritten in the form 

Wi} [ A j=l,2 J 

From Eq.(D we easily obtain 

9n{+) = e'-^^''^9r-li+)9r-li+) + 2e~''^gr-li+)9Vii-) + e'-'^:i:l(-)^7;i:l(-), 

9nH = e''-^^''^gr-i{+)9r-l{+) + 2e''' 9r-l{+)9r-l{-) + e-''^-''9r-l{-)9r-l{-)- 
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We introduce the following variable: 



For Xn we can then obtain the recursion relation 

= fix) = ^.^.(,-1) ^ + r 

where z = e^'^^, fi = e^^ and < a;„ < 1. The function f{x) is unimodal: it is continuous, 
continuously differentiable, and has one maximum x* in [0,1]. Note that f{x*) = 1 for 
any 7, h and T. This function is nonhyperbolic (hyperbolicity for ID maps means that 
1 < I/' I < 00 in all points) and maps the interval [0, 1] onto [z, 1]. 

Through x„, obtained by Eq.(|^), one can express the magnetization of the central base 
site: 

_ / s e^9:i{+) - e-'^g^iH _ e^^xl - 1 
" " e + e-N^(-) " e^'^xl + 1 ' 

and other thermodynamic parameters, so we can say that the x„ determines the states of 
the system. For example, at high values of temperatures (see next subsection) the recur- 
sion Eq.(^ tends to a fixed point and therefore the system has an appointed magnetization 
m. For the free energy function (using Eqs.(P|), (11), (1^), (0)), we obtain 



F„ = --J:i--lnt+'^ — ^^ln{zt^+2t+z){t^+2zt+l)-'^ — "^Inizt^+^t^+^zt+l), (8) 
3 2 2 3 

where t = fj.x^~^\ In deriving Eq.(P) we have used the following relation: 

F = Fcayley, (9) 

where Fcayiey is the free energy function of the triangular Cayley tree. It is easily seen 
that the expression for magnetization given by Eq. ( ^ , can be obtained by differentiating 
the free energy function of Eq.(P) with respect to the magnetic field h. It is interesting 



to mention that Eq.(^ is the generalization of the results obtained in Ref . |ll5|, |16 

3.THERM0DYNAMIC LIMIT 



Let us consider the magnetization of the central base site. In order to achieve ther- 
modynamic limit we tend the number of generation to infinity {n —>■ oo). If we set 7 = 4 
in Eqs.(^,(|^) and vary the temperature and magnetic field , then for sufficiently large T 
we see that function f{x) has a stable fixed point at every value of h and therefore map 
Xn = f{xn-i) attracts every point to x* = f{x*) in the thermodynamic limit [n —>■ 00). 
At such temperature the magnetization m, is a well defined function of h (see Fig. 2a) in 
thermodynamic limit. Then we lower T, at some values of h the point x* = f{x*) becomes 
unstable , but now two new stabile points xi,X2 arise in function f{f{x)) , which is an 
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attracting periodic orbit of period 2 for map x„ = f{xn-i) and we find that there is a 
single bubble in the plot of m versus h (Fig. 2b). As we continue to lower T, new bubbles 
are formed as parts of the old bubbles (Fig. 2c, attracting periodic orbits of period 2^^), 
and for still lower T s we reach a region where for intermediate values of h we have chaos, 
period-three windows, etc. (see Fig. 2d). 

One can say that the reason for this is that we have the presence of frustration effects 
and taking into account antiferromagnetic nature of the three-site interaction identifies two 
different magnetization in fixed field and temperature with magnetization of sublattice, as 



in antiferromagnetic phase and chaotic region with spin glass phase [^]. But as one can 
note from our calculation different magnetization corresponds to different sample with 
different number of generation. Therefore the appearance of a bubble dues to fact that 
our system does not have trivial thermodynamic limit and in that limit different samples 
with different number of generation are not macroscopically equivalent systems. 

Obviously magnetization can't be an order parameter in chaotic region. How we can 
identify the appearance of chaotic region? From the theory of dynamical system we know 
that Lyapunov exponents are good order parameter and become positive when attractors 
of the system become chaotic or strange]^. A statistical-thermodynamic method for 
the description of fluctuation of the Lyapunov exponents has been introduced, using a 
partition function[|l], §]. Later on we use thermodynamic formalism in order to characterize 
TSAI system in chaotic region. 

It is interesting to note that similar chaotic behavior has been found in other frustrated 
hierarchical lattices p!7|, 0. The Husimi tree like other hierarchical lattices are effectively 



infinite dimensional. But in Husimi or Bathe like lattices the number of neighborhood 
remain constant in contradiction to other hierarchical lattices for which in thermodynamic 
limit surface sites interacts with infinitely many neighborhood. Another difference is in 
nonequivalentce of the lattice sites. In Husimi or Bethe like lattices the sites lying deep 
inside lattices are equivalent []I9 . 



4.TSAI SYSTEM AND "THERMODYNAMIC FORMALISM" 

How to apply the "thermodynamic formalism" to the TSAI system? For this we need 
a natural partition and that is provided by the cylinders (we follow here Ref. |0]). 
The recursion function Eq.(|0) with 7 = 4 has the following form 



Simultaneously considering the following system of equation. 



l/(l)=Xo ' 



we obtain 

_ ^-2/3 ^ ^-8/3 _ 2^-1 
^~ 1 + - 2/X-5/3 ■ 



Xo = /i"'/'. (12) 
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For a crisis map (Eqs.(0), (|T2|)) we want to describe the scaling properties of the 
attracting set which in this case is the entire interval J: [xq, 1] (Fig.3). For an index n, 
/ is partitioned into 2" intervals or n-cylinders, these being the segments with identical 
symbolic-dynamics sequences of length n taken with respect to the maximum point x* = c. 
The inverse function of Eq.(0), h = f~^, has two branches, /i-i and hi as shown in (Fig.3) 
and the n-cylinders are all the nth-order preimages of /. The length of the cylinders is 
denoted by lei,e2,...t„ = h^^ o h^^ o . . . o h^^{I) where e G { — 1, 1}. 

The partition function Z{I3) is defined [H as 

Zn{(3)= E C....„= E e-'^'-'^i — , (13) 

where (3 G (— oo, oo) is a free parameter, the inverse "temperature". In the limit n ^ oo 
the sum behaves as 

^n(/5)=e-'^^^(^), (14) 
which defines the free energy, F{(3). The entropy S{X) is the Legendre transform 

5(A) = -/5F(/3) + A/3, (15) 

where the relation between A and /3 is found from 

A = ^(W/3)), m = S\X), (16) 

and these have the following meaning: In the limit n — > oo, e"'^'-'^^ is the number of 
cylinders with length / = e~"'^ or, equivalently, with Lyapunov exponent A, which we 
consider as an order parameter for TSAI model . The Hausdorf dimension of the set of 
points in / having Lyapunov exponent A is S'(A)/A 

By using the Eqs.(|10|), (p!2D, (p!3|), (|T^) we can numerically calculate the free energy at 



different value /i which is shown in Fig. 4. One can see from Fig.4a that the free energy 
have a singularity around /5c = —1 at low values of /i while at high values of the free 
energy have a singularity around jSc = 1 (Fig 4b), which shows the existence of the phase 
transition of first order in this regions of (3. 

How can the transition be determined accurately? 

Let us consider the characteristic Lyapunov exponent as an order parameter which 
will differ in the two phases. Fig. 5 shows this order parameter for different /i and sizes 
of the system, corresponding to n = 7,9,11,13. The curves converge towards a line 
and the result is the first order transition, whereas in a class of maps close to a; — ^ 
4a;(l — x) a phase transitions of the first order occur only with negative Pc- We would 
like to mention that it is hard to determine the critical "temperature" of phase transition 
by numerical methods with high precision. Consequently, the obtained value of critical 
"temperature" is approximate. In connection with this it is important to note, that there 
is a literature on phase transition in fully developed chaotic maps with neutral points 
leading to intermittent dynamics [Q, that claims the existence of a phase transition in a 
free energy at "temperature" = 1. 
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Large deviations of fluctuations of a Lyapunov exponents can be described by the 
entropy. To consider the above resuhs in terms of the entropy function S{X), let us first 
discuss the general appearance of the entropy function. First of all, it should be positive 
on some interval [Xmin, ^max]- The value A = ln2 must belong to that interval, which 
follows from the fact that the sum of the lengths of all cylinders on a given level is 1. 
Secondly it is often found that the values of Xmin and Xmax are given by the logarithms 
of the slopes at the origin . 

The precise form of the entropy function is, as mentioned in the introduction, not easy 
to obtain with great accuracy. The existence of a first order phase transition implies that 
there should be a straight line segment in S{X) and the slope of the line equal to Pc- This 
scenario is seen in Fig. 6. The curve in the figure corresponds to n = 13 Of course, with 
the finite-size data, it is impossible to determine the straight line segment in ^(A) and 
with increasing n, the straight line will be increased. 



6. CONCLUSION 



In this paper we have investigated the TSAI model by approximating it with a Husimi 
tree structures in an external magnetic field, and a strong connection with results from 
the theory of dynamical systems including chaos has been pointed out. 

An exact connection between that statistical system and fully developed chaotic at- 
tractors is obtained. Is is shown that in the chaotic region the magnetization, which are 
the order parameter for ferromagnetic phase, must be replaced by Lyapunov exponent in 
order to characterize TSAI model in chaotic region. The chaotic properties of the anti- 
ferromagnetic multisite system is described via the invariants characterizing a fractal set 
(e.g. a strange attractor). It is shown that this system in the chaotic region displays a 
phase transition. It is, in general, hard to determine such behavior unambiguously by nu- 
merical methods if one does not have further arguments or exact solutions. In this paper 
we have considered the characteristic Lyapunov exponents as an order parameter, which 
will differ in two phases. We've managed to calculate this order parameter as a function 
of the "temperature" and approximately showed where the phase transition will occur, 
since it describes transitions in the distribution of the characteristic Lyapunov exponents. 
This phase transition in terms of free energy and entropy function is also analyzed. 

It is interesting to note that for 7 = 3 the above situation changes dramatically and 
it will be very difficult to study numerically the chaotic region of TSAI system in terms 



of the "thermodynamic formalism" [12 



On other hand, the study of chaotic statistical physical system has opened new chal- 
lenges for theories of stochastic processes, especially in the direction of stochasticity of the 
vacuum in QCD [^, In this direction interesting results for the Z{Q) gauge model 



with a double plaquette representation of action on the fiat and generalized Bethe lattices 
were obtained[^, Note that it is possible to get a chaotic region in the Z{Q) gauge 
model with a three plaquette representation of the action. The detailed investigations of 
these questions will be published elsewhere. 
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Figure Captions 



Fig.l. Husimi tree with 7 = 4. 
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Fig. 2. Plots of m versus h' for different temperatures T (7 = 4). a - T = 3, b - 
T = 1.3, c - T = 1.15, d - T = 0.7. 
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Fig. 3. The function of Eq.(|l0D for values of /i = 105 and z given by Eq.(|T2]). 
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Fig. 4. The free energy F{(3), for different values of /x. a - = 5, b - = 105. 
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Fig. 5. The order parameter A(/3) calculated for different sizes of the system, corre- 
sponding to n — 9, 11, 13 and different /i. a-/x = 5, b-/i = 105. 
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Fig. 6. The entropy function corresponding to n = 13 for different /x. a- = 5, b - 
// = 105 . 
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